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Abstract The radical coupling reaction leading to the fi-O' 
quinonemethide intermediate of p-coumaryl alcohol was 
analyzed by semi-empirical molecular orbital calculation 
with MOPAC2002. By analyzing the radical monomer in a 
one-electron oxidation, the spin density of the unpaired 
electron at the 4-oxygen was less than half of the values at 
the C,, C 3 , C 5 , and Cp positions. By analyzing the transition 
state during the radical coupling reaction, the activation 
enthalpy was evaluated as 9.76kcal/mol, which corresponds 
to the activation energies for the propagation of common 
vinyl polymers. From the analysis of atomic interactions in 
the transition state, it was found that the activation enthalpy 
was largely composed of a high coulombic repulsion 
between Cp of the first monomer and the phenolic oxygen 
of the second monomer. After passing the transition state, 
the two radical monomers formed a metastable quinone¬ 
methide intermediate. The optimum conformation of the 
quinonemethide intermediate was formed from the meta¬ 
stable conformation through a second transition state with a 
small energy barrier. 

Key words Lignin • Radical coupling reaction ■ Transition 
state • Molecular orbital calculation ■ Activation enthalpy 


Introduction 

Lignin is formed by dehydrogenative polymerization 
(DHP) of monolignols oxidized enzymatically by the 
peroxidase/H 2 0 2 system. The overall reaction of DHP 
progresses in six stages: 
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1. Diffusion and collision of the monolignol and enzyme in 
the polysaccharide matrix. 

2. Penetration of monolignol into the active site of the 
enzyme. 

3. Formation of the radical molecules by one-electron 
oxidation of the monolignols, and their release from 
the enzyme. 

4. Formation of quinonemethide intermediates (generally 
fi-O', and /3-5' linkages) with coupling reaction of the 
two radical monomers. 

5. Stabilization of the a-carbon in the quinonemethides by 
the addition of a hydroxyl group for fi-O', resinol struc¬ 
ture for /3-/T, and coumaran structure for /1-5'. 

6. Further polymerization by repeating the first five 
stages. 

In recent years, quite a number of reports on the com¬ 
puter simulations of lignin formation have been published. 
For the first stage of DHP, Houtman and Atalla 1 simulated 
the diffusion of monolignols near the cellulose surface, after 
which Houtman 2 simulated the collision of two monolignol 
molecules by molecular dynamics calculations. For the 
third stage of DHP, Kobayashi et al. 3 simulated the activity 
of peroxidase to monolignols and related substances by 
semi-empirical molecular orbital calculations. The fourth 
stage was addressed by Elder and Ede 4 who simulated the 
theoretically feasible linkages of radical coupled intermedi¬ 
ates, and compared their reactivities by the heat of forma¬ 
tion. For the sixth stage, Glasser and Glasser 5 simulated 100 
units of lignin based on a statistical distribution of linkage 
types. However, the second and fifth stages of DHP have 
not been simulated yet. In other studies, a number of 
reports have emerged regarding lignin formation by ex¬ 
amining optimum conformations of final products after 
stabilization (fifth stage) by molecular orbital calculation. 
These reports were simulated with molecular mechanics, 6 
semi-empirical molecular orbital calculations, 7 ^ 9 or ab initio 
molecular orbital calculations. 10,11 However, these studies 
described the reactivity of DHP from the simulation of the 
stability of the final products and did not consider the reac¬ 
tion pathway. 
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Polymerization of lignin progresses through a transition 
state during the approach of two radical monomers (fourth 
stage), and the final products are formed following the for¬ 
mation of intermediates (fifth stage). We have focused on 
the former reaction, which determines the linkage types 
between two monomers. To search for a feasible pathway 
leading to coupled intermediates, the reverse reaction from 
coupled to uncoupled monolignols was computed. To sim¬ 
plify the calculation, p-coumaryl alcohol was selected as the 
model substance for monolignol, although it is an oversim¬ 
plification of lignin structure. The reason for the selection 
was the lack of conformationally complicated methoxyl 
groups , 12,13 which are present in coniferyl and sinapyl 
alcohols (constituents of common lignin). Furthermore, 
modeling of the coupling reaction was limited to formation 
of the fi-O' linkage, which is the main linkage type of native 
lignin. All modeling was conducted by semi-empirical mo¬ 
lecular orbital calculations rather than by ab initio methods 
due to the less demanding computer requirements for 
semiempirical methods. 


Methods 

Computation was carried out using WinMOPAC Ver. 
3.5.1a (distributed by Fujitsu). This software combined a 
graphical user interface for the development of the Z- 
matrix and MOPAC2002 Ver. 1.00 as the calculating 
engine . 14 The eigenvector following (EF) routine was used 
as the optimizer, and the gradient normal (GNORM) was 
set to 0.01. The Parametric Method 5 (PM5) was used to 
calculate the heat of formation ( H t ). The unrestricted 
Flartree-Fock Hamiltonian approximation (UHF) was used 
for the monolignol radical and the coupling reaction due to 
the presence of an unpaired electron. For other parameters, 
default values programmed in MOPAC2002 were used. 
The calculations proceeded as follows: 

1. Full optimization of p-coumaryl alcohol radical by taking 
into consideration the geometry of the alcohol hydroxyl 
group, which rotates freely. All bond lengths, angles, and 
torsions were optimized in the computation. 

2. Confirmation of the reactive atoms in the radical 
molecule by the localization of the unpaired electron 
estimated by the calculated spin density. 

3. A search for the conformers of fi-O '-linked 
quinonemethide intermediate by considering the geom¬ 
etries of C,-C„, C a -Cp, C /r C) 4 ., O 4 -C 4 ., and C r -C (/ bonds, 
and two alcohol hydroxyl groups. 

4. A search for the transition state candidates in the 
coupling reaction by disconnecting the fi-O' bond in the 
conformers of the quinonemethide intermediate. All 
bond lengths, angles, and torsions except the discon¬ 
nected bond length were optimized in the computation. 
The candidates were determined from the conformation 
giving maximal H f in each diagram of atomic distance 
and II,. 

5. Full optimization of the transition states (TS) from the 
Cartesian coordinate geometries of candidates as an ini¬ 




Fig. 1 . Carbon numbers in p-coumaryl alcohol radical (left), and 
quinonemethide intermediate of (S-)/f-0'-linked p-coumaryl alcohols 
(right) 

tial structure. Confirmation of obtained TS structures 
was done by a force calculation (FORCE) and an intrin¬ 
sic reaction coordinates calculation (IRC). The global 
TS structure in fi-O' coupling was determined as the one 
having the smallest energy of all TS structures. 

6 . Analysis of atomic interaction at the transition state by 
energy partitioning (ENPART). 

The carbon numbers in the molecules appearing in this 
report are defined as shown in Fig. 1. The definition 
of torsion angles, i.e., the relative positions of atoms, is 
presented in the text. 

The probability of minimum energy conformations was 
calculated by Boltzman’s distribution. The relative 
probability (p) was obtained from Eq. 1: 

p = e -[H t -«?)/RT (!) 

where R is the gas constant (1.98calK” 1 mol -1 ), T is the 
absolute temperature (298K in this report), and //[’ is the 
lowest H t value. 


Results and discussion 

Optimum conformation and reactive positions of p- 
coumaryl alcohol radical 

For optimum conformations of p-coumaryl alcohol radical, 
four minimum energy conformers were obtained. The opti¬ 
mum values of torsions and the geometries are shown in 
Table 1 and Fig. 2. Conformers a' and b' are symmetrical to 
a and b, respectively. The geometrical difference appeared 
in only the alcohol hydroxyl group, and the other geom¬ 
etries were almost identical. The energy difference between 
a and b was small and corresponded to 1:0.385 of the exist¬ 
ing probability. The lowest H { for the p-coumaryl alcohol 
radical was -36.76kcal/mol. Therefore, H t of the uncoupled 
state, i.e., two infinitely separated molecules, was evaluated 
as -73.52kcal/mol (twice the value for a single molecule). 
This corresponds to the initial state of the radical coupling 
reaction. 

The reactive positions of the p-coumaryl alcohol radical 
were analyzed from the spin population and atomic charge. 
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Table 1 . Optimum conformations of p-coumaryl alcohol radical 


Conformation 

Oy-CfCp-Ca (degree) 

H-O r -C,-Cp (degree) 

Hf (kcal/mol) 

P 

a 

24.6 

60.6 

-36.76 

1 

a' 

-24.6 

-60.6 

-36.76 

1 

b 

156.6 

-62.9 

-36.20 

0.385 

b' 

-156.6 

62.9 

-36.20 

0.385 


b 

Fig. 2. Geometries of p-coumaryl alcohol radical optimized by UHF/ 
PM5 protocol. Conformations correspond to those in Table 1. Confor¬ 
mations a' and b' are symmetries of a and b, respectively 


Table 2 shows the spin populations of a-spin and /1-spin, the 
electron density, the spin density, and the Mulliken charge 
for conformer a of the p-coumaryl alcohol radical. The spin 
density was found to be high at C 4 , C 3 , C 5 , Cp, and 0 4 
suggesting the localization of an unpaired electron at these 
atoms, and implying the occurrence of radical coupling 
reaction. The atoms having high spin densities were similar 
to the coniferyl alcohol radical. 913 Additionally, the re¬ 
active atoms were negatively charged, suggesting that the 
radical coupling reaction proceeds despite the coulombic 
repulsion. 

Figure 3 shows the distributions of the unpaired electron 
in p-coumaryl alcohol radical and the relative ratios pre¬ 
dicted from the spin densities at C q , C 3 , C 5 , Cp, and 0 4 . By 
assuming the coupling reaction occurs in accordance with 
these ratios, the radical coupling appears to occur predomi¬ 
nantly at C 4 , which has the highest spin density. However, 
coupling at C, would be difficult due to steric hindrance. 
Other couplings, fi-fi' and /J-5' may easily occur due to the 
high ratios at C 5 and Cp. The formation of fi-O' coupling 
may be comparatively difficult due to the low ratio at 0 4 . 

Unfortunately, semi-empirical molecular orbital calcula¬ 
tion for the unpaired electron is less reliable than advanced 
methods including electron correlation theory, e.g., the den¬ 
sity functional models, the Mpller-Plesset models, and more 
advanced methods. The details of comparisons between dif¬ 
ferent monolignols and between different couplings will be 
discussed in subsequent reports. 

Conformations of quinonemethide intermediate 

In the fi-O' quinonemethide intermediate, many conforma¬ 
tions are possible due to the rotation at the C,-C„, C a -Cp, C p- 
C) 4 .. 0 4 -C 4 ., and C r -C (/ bonds, and the two alcohol hydroxyl 
groups. Thus, the following steps were taken to arrive at the 
minimum energy conformations. 5- and /Cform enanti¬ 


omers at the chiral center Cp are present in the fi-O '-coupled 
intermediate. But as MOP AC cannot analyze the difference 
between the enantiomers, only the 5-form was analyzed in 
this study. For the Q-C,* bond, the position of Cp ((j)p is 
considered to be mostly in the same plane as the first 
monomer’s aromatic ring due to the presence of a double 
bond at Cf-C a . For the C a -Cp bond, the position of C r 
(QrP-a- 1 ) ' s considered to be anti to C, due to the monomer 
structure. Simultaneously, the position of C) 4 . is restricted 
by the position of C 7 due to the sp 3 orbital of Cp. For the 
Cp-C) 4 . bond, the position of C 4 < can be defined as the torsion 
angle against C a (</> 4 -. 0 ’.p-d, and would be restricted at three 
gauche positions against C a , C r , and Hp. For the 0 4 -C 4 < 
bond, the aromatic ring of the second monomer is mostly 
restricted to be parallel to the aromatic ring (first 
monomer)-C a -Cp plane by steric hindrance. For C r -C rf , Cp 
is limited to two positions against the aromatic plane of the 
second monomer (definable as </v 4 '.ovi) analogous with the 
conformation of the monomer. In addition, four geometries 
of the alcohol hydroxyl group (a, a', b, and b') shown in the 
monomer structure must be considered for each monomer. 
These 96 possible conformers (Ixlx3xlx2x4x4) were 
used as the initial structure for full optimization. In addi¬ 
tion, the feasible transition states were derived by discon¬ 
necting the fi-O' bond from all optimum conformers of 
quinonemethide intermediate. Table 3 shows the low H f 
conformations of the quinonemethide intermediate of 
(5-)/3-0'-linked p-coumaryl alcohols, as well as the low en¬ 
ergy transition states. Typical examples of the geometries, 
QM1 and TS1, are shown in Fig. 4. 

In low energy conformers of the quinonemethide inter¬ 
mediate (QM1-8), the positions of C 4 . were similar and anti 
to Hp in the 0 4 -C r , bond. Furthermore, rotation of the 
second monomer’s aromatic ring at the C 4 -0 4 . bond was 
restricted within the range of 70°-90°; the absolute value in 
C 3 < shows the torsion of the two aromatic planes, and the 
sign shows the direction of Cp. The y-hydroxymethyl group 
was restricted by a and b forms but not a' and b' (see 0 7 H r ). 
These conformational restrictions could possibly be due to 
the united effects of hydrogen bonding and steric hindrance 
between the hydroxyl group and the aromatic ring. The 
influence of the Cp position on H t was minimal because of 
the appearance of both positive and negative signs of C 3 . in 
Table 3. Furthermore, the y'-hydroxymethyl group was re¬ 
stricted to a and a' forms (see O r TIy), which are optimum 
conformations that are analogous to the results of the 
radical monomers shown in Table 1. These results suggest 
negligible atomic interactions with the first monomer’s 
atoms due to their distance of separation. 






Table 2. Spin population and charge density of each atom in conformer a of p-coumaryl alcohol 
radical computed by UHF/PM5 protocol 
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Atom Spin population Mulliken atomic charge 



a-spin 

/1-spin 

Electron density 

Spin density 


Q 

2.342 

1.673 

4.015 

0.669 

-0.015 

C 2 

1.825 

2.294 

4.119 

-0.469 

-0.119 

C 3 

2.398 

1.849 

4.247 

0.549 

-0.247 

C 4 

1.741 

1.912 

3.653 

-0.170 

0.347 

C 5 

2.394 

1.854 

4.247 

0.540 

-0.247 

C 6 

1.826 

2.293 

4.119 

-0.467 

-0.119 

C„ 

1.824 

2.309 

4.132 

-0.485 

-0.132 

C, 

2.417 

1.826 

4.244 

0.591 

-0.244 

c r 

1.978 

2.033 

4.011 

-0.054 

-0.011 

o 7 

3.214 

3.210 

6.424 

0.004 

-0.424 

o 4 

3.323 

3.048 

6.371 

0.274 

-0.371 

H 2 

0.431 

0.404 

0.835 

0.028 

0.165 

h 3 

0.386 

0.420 

0.805 

-0.034 

0.195 

h 5 

0.385 

0.419 

0.804 

-0.034 

0.196 

h 6 

0.429 

0.402 

0.830 

0.027 

0.170 

H„ 

0.429 

0.400 

0.829 

0.029 

0.171 


0.400 

0.435 

0.834 

-0.035 

0.166 

H r 

0.436 

0.413 

0.849 

0.023 

0.151 

h 7 

0.456 

0.443 

0.898 

0.013 

0.102 

H 0 

0.365 

0.365 

0.731 

0.000 

0.269 

Total 

29.000 

28.000 

57.000 

1.000 

0.000 


Electron density: sum of a-spin and /1-spin.; spin density: excess of a-spin over /1-spin 



0.410 0.821 0808 1 0.884 


Fig. 3. Distribution of unpaired electron in p-coumaryl alcohol radical. 
The values indicate the relative ratios predicted from the spin densities 
at Q, C 3 , C 5 , Cp, and 0 4 


Conformations and intermolecular interactions between 
radical monomers at the transition state 

In the low energy conformations of the transition state 
(TS1-8) shown in Table 3, the positions of 0 4 - are similar; 
the C^-O^ distances were about 1.93 A, and the O^-C^-C,* 
bond angles were about 100° for all transition states. These 
similarities would be governed quantum mechanically by 
the sp 3 orbital of Cp and 0 4 .. In these transition states, the 
positions of C 4 - of the second monomer were similar, and 
restricted to be anti to C r The structural difference between 
these transition states and coupled intermediates may be 
due to the strong hydrogen bonding of the y-hydroxyl group 
to 0 4 - of the first monomer, and the small interaction with n- 
electrons in the second monomer’s incomplete aromatic 
ring. 

For the lowest energy conformation (TS1), only 
one imaginary vibrational frequency was obtained at 
-649.88 cm -1 in the direction of approach of Cp and () 4 .. This 
conformation was verified as the transition state because 


Table 3. The major conformers extracted from 96 initial structures of 
quinonemethide intermediate of (5-)/TO' linked p-coumaryl alcohols 
(QM1-8) and transition states toward quinonemethide intermediate 

(TS1-8) 


No. 

Conformation 




H l (kcal/mol) 

P 

c/ 

c 3 - b 

0,H, C 

O r .H/ 


QM1 

H, 

73.3 

a 

a' 

1.414 

-89.33 

1 

QM2 

H, 

73.5 

a 

a 

1.414 

-89.33 

0.993 

QM3 


-112.0 

a 

a' 

1.414 

-89.27 

0.905 

QM4 

H„ 

-112.1 

a 

a 

1.414 

-89.27 

0.901 

QM5 


88.9 

b 

a 

1.416 

-88.89 

0.479 

QM6 


89.3 

b 

a' 

1.416 

-88.88 

0.465 

QM7 

H, 

-97.0 

b 

a 

1.416 

-88.85 

0.442 

QM8 

H, 

-96.8 

b 

a' 

1.416 

-88.84 

0.438 

TS1 

c r 

-88.1 

a 

a' 

1.931 

-63.76 

1 

TS2 

C r 

-90.3 

a 

a 

1.935 

-63.76 

0.995 

TS3 

C r 

92.6 

a 

a 

1.930 

-63.63 

0.800 

TS4 

c. 

-93.1 

b 

a 

1.936 

-63.40 

0.537 

TS5 

C r 

-93.0 

b 

a' 

1.933 

-63.34 

0.491 

TS6 

c r 

-87.6 

a 

b' 

1.931 

-63.33 

0.484 

TS7 

C r 

89.9 

a 

a' 

1.934 

-63.32 

0.473 

TS8 

C? 

87.0 

b 

a 

1.931 

-63.21 

0.391 


“Anti -atom of C 4 . with Cy0 4 . bond 
b Torsion angle of C 3 .-C 4 .-0 4 .-C^ 

‘Corresponding to the conformation name of monomer shown in 
Table 1 

‘‘Distance between Cp and 0 4 . in angstroms 


only one imaginary frequency was obtained. Also, with the 
direction indicating the approach of Cp and C) 4 ., this transi¬ 
tion state was verified to be reactive toward the j3-0' 
linkage. The activation enthalpy (//,,) was calculated 
as 9.76kcal/mol by using the equation // a = // f (TSl) - 73.52. 
This value corresponds to the activation energies of radi¬ 
cal polymerization of common vinyl polymers, e.g., 










132 




Fig. 4. Typical examples for geometries of the quinonemethide inter¬ 
mediate of (S-)/5-0'-linked p-coumaryl alcohols (left, QM1) and the 
transition states during radical coupling reaction (right, TS1) 


Table 4. The major intermolecular interactions between radical 
monomers at TS1 


Atoms 


J (eV) 

K (eV) 

C(eV) 

J + K+C (eV) 

c p 

o 4 . 

-2.7221 

-0.9262 

1.4445 

-2.2038 

H 0 

O, 

-0.0424 

-0.0045 

-0.5444 

-0.5914 

C 4 

Oy 

0.0000 

0.0000 

-0.3496 

-0.3496 

C 4 

o 4 . 

0.0000 

-0.0006 

-0.3344 

-0.3351 

o 4 

H 0 - 

0.0000 

0.0000 

-0.3198 

-0.3198 

H„ 

O, 

-0.0058 

-0.0010 

-0.2905 

-0.2973 

O r 

C 4 - 

-0.0004 

-0.0011 

-0.2891 

-0.2906 

c 3 

Oy 

0.0000 

0.0000 

0.2509 

0.2509 

o 4 

c„ 

0.0000 

-0.0001 

0.2574 

0.2574 

C 4 

Ho- 

0.0000 

0.0000 

0.2655 

0.2655 

C 3 

o. 

-0.0001 

-0.0016 

0.2674 

0.2657 

o 4 

O, 

0.0000 

-0.0008 

0.2895 

0.2887 

c„ 

o 4 - 

0.0301 

-0.0461 

0.3894 

0.3734 

o 4 

Oy 

0.0000 

0.0000 

0.4034 

0.4034 

Or 

O, 

0.0044 

-0.0033 

0.6922 

0.6933 


Negative and positive values mean attractive and repulsive inter¬ 
actions, respectively 

J, Resonance energy; K, exchange energy; C, coulombic interaction 


7-8kcal/mol for the propagation of polystyrene, 16 and is 
much lower than the activation energies of the decomposi¬ 
tion of common initiators, e.g., 26-34kcal/mol for 2,2'- 
azo-bis-isobutyronitrile and benzoyl peroxide. The entropy 
term may be assumed to be negligible, in which case the 
radical coupling reaction of lignin will be as fast as the 
propagation of vinyl polymers. 

The major intermolecular atomic interactions in TS1 
were calculated and the results are shown in Table 4. The 
largest total interaction was, of course, that at C^-O^ with 
-2.2038 eV of attraction. This atomic interaction contained 
large coulombic repulsion component caused by the nega¬ 
tive charges of both C /( and C 4 -. The excess attraction of 
resonance and exchange energies were caused by the ap¬ 
proach of the two atoms having high spin densities. For 
other atomic pairs, the total interactions were lower than 
1/3 of that of C«-0 4 ., and mostly due to coulombic interac¬ 
tion including hydrogen bonding. Additionally, this analysis 
revealed that the alcohol hydroxyl group contributes con¬ 
siderably to the optimum conformation of the transition 
state. In other words, the kind of substituent at C,,of mono¬ 
mers may affect the coupling reaction, e.g., carboxyl or 
aldehyde groups. 



Fig. 5. The heat of formation (filled squares) and the spin densities at 
Co (open circles) and 0 4 . (filled circles) during radical coupling reaction 
of p-coumaryl alcohols as a function of the Cp-C 4 < distance 


Figure 5 shows the changes of H t and spin densities at Co 
and C 4 - as functions of the atomic distance between them; 
this diagram was drawn by replotting the geometrical and 
orbital data on the reaction coordinate computed by IRC 
for the structure of TS1. In the coupling reaction, //, values 
increased as the two monomers approached each other, and 
reached a maximum at the transition state. During the ap¬ 
proach of the two monomers, both spin densities increased, 
especially at 0 4 -, indicating that the localization of both 
unpaired electrons changed due to the resonance of C p and 
0 4 -. The variation of spin densities during the reaction 
means that the radical coupling reaction cannot be ex¬ 
pressed by only the spin density in the isolated monomer 
given in Table 2. After passing the transition state, H t de¬ 
creased and reached the value of the quinonemethide 
coupled form (QM TS1 ), -89.33 kcal/mol. Because //, of 
QM TS1 is lower than that of the uncoupled state, it is obvious 
that this reaction can proceed. Simultaneously, the spin 
densities suddenly decreased and finally disappeared in the 
coupled state. Consequently, it was found that the 
two monomers formed the coupled state, QM TS1 , by IRC 
analysis. 


Transition from QM TS1 to QM1 

The coupled state, QM TS1 , is metastable and different from 
the lowest energy conformer of quinonemethide, QM1. 
QM TS1 and QM1 were structurally different at the C 4 . posi¬ 
tion, but not at C 3 .. Thus, the energy barrier in transition 
from QM TS1 to QM1 with rotation of the aromatic ring of 
the second monomer at the C^-O,,. bond was analyzed. Be¬ 
cause the analysis could not be achieved by UF1F protocol, 
it was performed by the restricted Hartree-Fock protocol in 
this study. Estimated values of H f were -87.27 kcal/mol 
(-88.01 kcal/mol by UFIF) for QM TS1 , -88.49 kcal/mol 
(-89.33 kcal/mol by UHF) for QM1, and -86.69 kcal/mol 
for the transition state between them. Therefore, // a 
for the transition from QM TS1 to QM1 was estimated as 
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Transition state (T S1) 

/ A • 

/ 

«■ 

_ 

Uncoupled state 

Transition state 


--- 

- 

Metastable state 

j Quinonemethide 

T ° intermediate (Q M 1) 


Radical coupling reaction —* 

Fig. 6. Overall energy change in radical coupling reaction of two p- 
coumaryl alcohols toward p-O' quinonemethide intermediate. The ver¬ 
tical axis is the energy difference from the heat of formation of two 
uncoupled monomers 


0.58kcal/mol. The extremely small energy barrier indicates 
that this transition easily progresses. 

The existence of different conformations between 
estimates from the transition state and the optimized 
quinonemethide is important for analysis of the overall 
reaction. This means that the thermodynamic and kinetic 
estimations may lead to different results. In other words, 
the prediction of the reactivity based on only the conforma¬ 
tion of products is insufficient. 


The overall radical coupling reaction of p-coumaryl 
alcohol leading to /j-O'-quininemethide intermediate 

Figure 6 shows the overall radical coupling reaction of p- 
coumaryl alcohol. The energy level based on the two un¬ 
coupled monomers (twice the value of conformer a) 
decreased to -15.81 kcal/mol at QM1. During this process, 
two transition states were detected by the molecular orbital 
calculation. One was the transition state in coupling the two 
monomers with // a = 9.76 kcal/mol. The other transition 
state was the conformation change from the metastable 
quinonemethide to the optimum state. However, the energy 
barrier of the second transition was extremely low. Thus, 
the rate-determining step in the radical coupling reaction 
was considered to be the first reaction, i.e., bond formation 
between C» and C 4 <. It is concluded that the radical coupling 
reaction occurs by kinetic control. 
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